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Characterization and Photometric Membership of the 
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ABSTRACT 

Open clusters belonging to star-forming complexes are the leftovers from the initial 
stellar generations. The study of these young systems provides constraints to models 
of star formation and evolution as well as to the properties of the Galactic disc. We 
aimed at investigating NGC 1981, a young open cluster in the Orion Nebula Region, 
using near-IR and BV{RI)c photometric data. We devised a method that accounts for 
the field contamination and allows to derive photometric membership for the cluster 
stars. A new cluster centre was determined by Gaussian fittings to the 2-D stellar 
distribution on the sky, and has been used used to obtain the radial stellar density 
profile and the structural parameters. Mass functions were computed for stars inside 
the cluster limiting radius and total mass estimated from them. Although more easily 
distinguished by its grouping of 6 relatively bright stars, an underlying population 
of faint pre-main sequence stars is evident in the cluster area. We showed that this 
population is related to the cluster itself rather than to the nearby Orion Nebula 
cluster. Additionally a fraction of the cluster low mass stars may have been evaporated 
from the region in its early evolution leading to the present sparse, loose structure. 
The estimated parameters of NGC 1981 are core radius Re = 0.09 ± 0.04 pc, limiting 
radius Rum = 1.21 ±0.11 pc, age t = 5±1 Myr, distance modulus {m — M)o = 7.9 ±0.1 
(380 ± 17 pc), reddening E{B - V)= 0.07 ± 0.03 and total mass m = 137 ± 14 M©. 

Key words: open clusters and associations: general - open clusters and associations: 
individual: NGC 1981 - galaxy: stellar content - stars: pre-main-sequence 



1 INTRODUCTION 

The open clusters' fundamental physical parameters are 
important pieces of information for studies o n the forma- 
tion a nd evolution of the G a lactic disc (e.^., Piatt i et al.l 
I1995I : IJacobson eraD I2OO8I : ISestito et all I2OO8I I and as 
groun ding tests for star formation and evolution mode ls 
fe.g.. iLandm et all 1 1996: DAntona '20021: ISiess e t"aL 2000). 
iDias et all (|2002l , hereafter DAML02) have summarized the 
available information on open clusters in a major catalogue 
which continuously expands. In version 2.1Clj of this cata- 
logue, 1787 open clusters are listed out of which 54 per cent 
have known distance, age and reddening, 24 per cent have 
also proper motion and radial velo city and only 10 p e r cent 
have metal abundance determined. iKharchenko et al.l (|2005l l 
have created an open cluster catalogue with uniformly deter- 
mined parameters, among them cluster radius, core radius 
and age, by means of an automated computational algorithm 



(The Catalogue of Open Cluster Data, COCD fl. This cata- 
logue has been used to investigate the local (< 1 kpc) pop- 
ulation of star cluster complexes revealing new constraints 
on th e Galaxy's structure and kinematics (|Piskunov et al.l 
I2OO6I ). 

Sky surveys like 2MASS (Skru tskie et all I2OO6I I pro- 
duced large amounts of near-IR data and have contributed 
to the discovery of even more objects. Additionally, sev- 
eral studies have been benefited from 2MASS database by 
employing near-IR photometric analyses of structural and 
populational properti e s of Milky Way star clusters (e.g., 
ISant os Jr. et all l2005l : IPavani fc Bical l2007l : iBonatto et aJ] 
2006). Because the 2MASS database covers the whole sky, 
it allows data extraction from spatially unlimited regions. 
Also, near-IR wavelengths are particularly sensitive to dis- 
criminate cluster stars fr om the contaminating field for 
young stellar systems fe.g. JSantos Jr. et al.ll2005l V 

NGC 1981 (OCl 525, C 0532-044) is a young, sparse clus- 
ter located ^1° North from the Orion Nebula at Galactic 
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coordinat e s I = 208.09° and h = -18.98°. As noted by 
ISharplessI ((1953), ^he Orion Nebula seems to he separated 
from NGC 1981 only as a result of a somewhat nearer ob- 
scuring cloud. NGC 1981 depicts the northern en d of Orion's 
Sword. According to lSubramaniam et al] (| 19951 ). NGC 1981 
and Cohinder 70 constitute a possible binary cluster be- 
ing less than 20 pc apart. The cluster belongs to Gould's 
Belt, a planar distribution of O and B stars inclined ~ 20° 
with respect to the Galactic plane (Poppel |200lL an d ref- 
erenc es therein). In a detailed study of Gould's Belt (|Leshl 
[1968'), the brightest stars in the cluster (HD 37016, HD 37017 
and HD 37040) had their spectral types determined (B2.5V, 
B1.5V, B2.5IV, respectively) and photometric distances es- 
timated (all at 494 pc). 

Subsequent works emphasized kinematical properties of 
NGC 1981 stars, among other clusters, aiming at obtaining 
dynam ical and structural characteristics of our Galaxy (e.g., 
lHronlll987l ). More recently. iKharchenko et al.1 (|2005l ) per- 
formed an analysis of the cluster as a stellar system deriving 
reddening E{B - V)= 0.05, distance modulus {V - M^) = 
8.16 {d = 400 pc), logt(yr) = 7.50, cluster radius R = 0.25° 
(1.7 pc) and core radius Re = 0.13° (0 .9 pc). We should 
note, however, that the age estimated bv IKharchenko et al.l 
(|2005[ ) relies upon a single star, as a consequence of the 
method applied, which progressively removes stars with low 
membership probabilities from the CMD by means of an 
iterative process. 

Most studies of the cluster are inserted in more gen- 
eral analyses of the Orion star formation complex. It is part 
of Orion OBI association, subgroup c, which is 2 to 6Myr 
old and located closer to us (^400 pc) th an the younger 
Orion Nebula by at least 10 pc (^Ballvl 120081 ). However, such 
a division in subgroups has be en questioned in fav our of 
a continuous star forming event rMuenc h et al.ll2008r ). Since 
NGC 1981 bright stars are early B spectral types, a plausible 
evolutionary sequence would entail supernovae explosions 
from 0-type progenitors causing the compression of the in- 
terstellar medium and formation o f the present youn ger pop- 
ulations. On the kinematical side, Dias et al.l (|200 6) derived 
membership probabilities on the basis of UCAC2 proper mo- 
tions aided by a statistical method for 158 stars within 15 
arc min of the c l uster centre . 

lElias et al.l (|2009|) investigated hierarchical star forma- 
tion in Gould's Belt based on the spatial and kinemati- 
cal distribution of star clusters, either bound star systems 
or transient stellar condensations with a mean lifetime of 
~10Myr or less. NGC 1981 was classified by them as a 
bo und cluster in view of it s age (~ 30 Myr) as determined 
bv IKharchenko et al.l (|2005l ). In the present work, we argue 
that NGC 1981 is indeed younger, as has been shown to be 
the case for Platais 6 and NGC 2546, which also belong to 
Gould's Belt (JEhas et al.ll2009l ). 

Despite the efforts to better understand the astro- 
physical processes occurring in the Orion Nebula and 
its surroundings and their connection with Gould's Belt, 
N GC 1981 have not been ^i ven particular attention (except 
for IKharchenko et al.l l2005l ) and, as a consequence, its pa- 
rameters can be improved. In this work, we present Johnson- 
Cousins BV{RI)c photometric observations carried out at 
the Observatorio do Pico dos Dias (Itajuba, Brazil) along 
with near-IR photometric data extracted from 2MASS sur- 
vey to study the open cluster structural properties and stel- 



lar content, aiming at a better determination of radius, age, 
distance, reddening and mass. 

This work is structured as follows. Sect. 2 describes the 
observational data employed. Sect. 3 presents the procedure 
to obtain the cluster centre, a fundamental step to build a 
reliable radial density profile, discussed in Sect. 4. Colour- 
magnitude diagrams are employed in Sect. 5, where a sta- 
tistical decontamination algorithm is applied to the cluster 
CMDs, and in Sect. 6, where a comparison of data with 
isochrones allows a determination of the cluster's astrophys- 
ical parameters. In Sect. 7 the radial mass function of the 
cluster is analysed and in Sect. 8 the stellar density charts 
used to map the spatial distribution of stars and their fre- 
quency in different CMD regions are described. Further dis- 
cussions and conclusions are given in Sect. 9. 



2 DATA 

2.1 Near- infrared 

Vizieiu was used to extract near-infrared photometric data 
from 2MASS in circular fields centred on equatorial coor- 
dinates (J2000) a = 05^35"'09' and 6 = -04°25'54'^ the 
cluster centre as taken from DAML02. The data encom- 
passes the point sources within i? < 55 arcmin, correspond- 
ing roughly to four times the estimated visual radius of the 
cluster. For decontamination purposes, a comparison field 
was extracted 1° northwest of the cluster, in a region with 
the same area and similar absorption as deduced from near- 
infrared (2MASS) and mid-infrared (IRAS) images. This 
choice was preferred over an annular field extraction cen- 
tred in the cluster because of the nebulosity associated with 
NGC 1977 towards the south. Figure [1] shows an image of 
NGC 1981 and the nearby region from IRAS 60 ^im band. 
The selected comparison field and other objects are also 
shown. 

As photometric quality constraint, the 2MASS ex- 
tracted data were restricted to unambiguous point sources 
brighter than specified on the 99 per cent Point Source Cat- 
alogue completeness limit (15.8, 15.1 and 14.3 mag for 10- cr 
at J, if, and Kg, respectively). 



2.2 Optical 

CCD BV{RI)c observations were carried out with the 0.6-m 
(lAG) telescope at the Observatorio do Pico dos Dias (Lab- 
oratorio Nacional de Astrofisica, Itajuba, Brazil) during the 
night of September 21, 2000. The CCD detector employed 
was a Ixlk SITe SI003AB (# 101) configured to yield read- 
out noise 5.5 e~ and gain 5.0e~ADU~^. The plate scale was 
1.22 arcsec. pixel" ^. The images were obtained using a focal 
reducer to cover a field-of-view of 21 x 21 arcmin and were 
taken with a mosaic setup on 6 adjacent positions around 
the cluster centre in order to cover the full extension of the 
object. To transform instrumental magnitudes into the stan- 
dard system, 15 objects from Landolt (1992) standard stars 
catalogue were also observed on the same night, with 2 stan- 
dard stars being used for atmospheric extinction correction. 
Details of the observations are shown in Table [U where the 
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Table 1. Observation log for NGC 1981 



Figure 1. NGC 1981 and the nearby region from IRAS 60 ixm 
band. Other objects in the area and the selected comparison field 
are also indicated. 



observed fields are listed along with their exposure times, 
typical seeing values and airmasses. 

The data were reduced using the standard iraf rou- 
tines for aperture photometry. Besides the standard reduc- 
tion steps, the 6 adjacent fields have been aligned and re- 
current stars were selected according to their distance to the 
image centre in each filter. Specifically, for every recurrent 
star, we kept only the measurement of the star closest to 
the image centre. This choice avoided the vignetting effect 
introduced on the images borders by the focal reducer. The 
calibration equations used to transform to the standard sys- 
tem were: 



h = B + bi + b2Xb + b3{B-V) + bA{B-V)Xi, (1) 

V = V + vi+V2X^+vs{B-V)+v^{B-V)X^ (2) 

r = R + ri+r2Xr + r3{V -R) + r^{V -R)Xr (3) 

i = I + ii + ^2X^ + isiV - I) + i^iV - I)X^ (4) 

where B, V, R and / are the standard magnitudes and 6, v, 
r and i are the instrumental magnitudes; Xn is the airmass 
on each filter. The subscript coefficients i — 1...4, were ob- 
tained from a 2-pass interactive fitting on each filter with the 
extinction coefficients (subscripts 2 and 4) being fitted first 
and then used as constants on a second fit of the zero point 
and colour terms (subscripts 1 and 3 respectively). The de- 
rived coefficients and the RMS deviation of the fittings to 
the catalogued magnitudes of the standard stars are shown 
in Table [21 

As photometric quality constraint the employed data 
set was that restricted to an error less than 0.04 mag, corre- 
sponding to the higher deviations found on the calibration 
fit. This value implies a magnitude limit of approximately 
15 mag and leads to an average magnitude uncertainty of 
0.05 mag (taking into account the calibration deviation). 

The following analysis involving spatial information 



Cluster/Region 


Filter 


Exposure (s) 


Seeing ('') 


Airmass 




B 


1 15 50 


3.4 


1.2 


NGC1981/CE 


V 


1 20 


2.5 


1.2 




R 


1 20 


2.4 


1.2 




I 


1 20 


3.9 


1.2 




B 


1 15 50 


3.1 


1.1 


NGC1981/CW 


V 


1 20 


2.8 


1.1 




R 


1 20 


2.7 


1.1 




I 


1 20 


3.4 


1.1 




B 


1 15 50 


2.7 


1.3 


NGC1981/NE 


V 


1 20 


2.6 


1.3 




R 


1 20 


2.3 


1.2 




I 


1 20 


3.8 


1.2 




B 


1 15 50 


3.1 


1.1 


NGC1981/N 


V 


1 20 


2.6 


1.1 




R 


1 20 


2.6 


1.1 




I 


1 20 


3.7 


1.1 




B 


1 15 50 


3.3 


1.1 


NGC 1981/NW 


V 


1 20 


2.9 


1.1 




R 


1 20 


2.8 


1.1 




I 


1 20 


4.1 


1.1 




B 


1 15 50 


3.3 


1.1 


NGC1981/S 


V 


1 20 


2.6 


1.1 




R 


1 20 


2.7 


1.1 




I 


1 20 


3.7 


1.1 



Table 2. Coefficients and RMS deviation of the calibration fit. 



V 



R 



1 


5.79±0.01 


5.35±0.01 


4.53±0.02 


5.05±0.01 
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0.28±0.06 


0.17±0.04 
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4 
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0.02 
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were developed using 2MASS data due to its large extrac- 
tion area, good astrometry and the poor seeing of the opti- 
cal data, which was employed mostly for CMD comparisons 
with the near-IR data. 



3 NEW CENTRE DETERMINATION 

The sparse aspect of NGC 1981 makes difficult to estimate 
its centre, which is an essential step before a reasonable ra- 
dial density profile can be obtained. Catalogued centre val- 
ues are mainly intended for identification purposes and are 
generally imprecise. Therefore, we have estimated the centre 
of NGC 1981 by first selecting a region (see Fig|2]) around the 
coordinates given by DAML02. The search area was limited 
to the sky around the catalogued centre and did not include 
stars southward of ^ = —4.55° to avoid the nebulosity and 
contamination from NGC 1977. 

Furthermore, the selected region was divided into bins 
of right ascension and declination and star counts were made 
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Figure 2. Sky chart of NGC 1981 showing the region used for 
centre determination (rectangle), the hterature centre (cross), the 
calculated centre (circle) and the data extraction area for the 
decontamination method (dotted circle). 
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Figure 4. Relation between determined coordinates and the bin 
size (top). Histograms showing the most recurrent right ascension 
and declination for the cluster centre (bottom). 



that the calculated centre is 5.32 arcmin north and 0.27 
arcmin east from the literature centre. The area used for 
centre determination and the data extraction region for the 
decontamination method (see Sect. 5) are also shown. 
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Figure 3. Gaussian fittings for the determination of central co- 
ordinates. Bin size on RA and DEC is 3.0 arcmin. Error bars 
correspond to Poisson statistical errors. 



inside them. We used these star counts to build spatial pro- 
files and fitted a Gaussian function to the stars distribution 
on RA and DEC, adopting the centre of the fitted Gaussian 
as the cluster centre coordinate. Fig. [3] shows the Gaussian 
fitting obtained for bin size 3.0 arcmin and the resulting 
coordinates of the centre. 

This procedure was applied to different bin sizes ranging 
from 0.25 — 4.00 arcmin and the resulting centre coordinates 
were used to create histograms showing the most recurrent 
RA and DEC values. We adopted these values as the coordi- 
nates of the cluster centre. We also found no trends between 
the coordinates and the bin size. Figure 2] shows this rela- 
tionship and the histograms defining the cluster centre. 

The adopted coordinates of the cluster centre are 
a = 83°46'59'' = 5^35"^08' and 6 = -04°20'35''. The 
histograms bin width were chosen based on the standard 
deviation of the centre coordinates and accounts for an un- 
certainty of approximately 5 arcsec on a and S. Fig.[2]shows 



4 RADIAL DENSITY PROFILE 

The structural properties of star clusters can be derived by 
means of a projected density profile such as the radial den- 
sity profile (RDP). The RDPs of star clusters usually follow 
an analytical profile and can be described by functions with 
different parameters related to the cluster structure. These 
parameters include the tidal radius, core radius, core stel- 
lar density and background star density. However, young 
and sparse open clusters like NGC 1981 usually show little 
overdensity in comparison to the field making difficult the 
fitting of analytical functions and so the determination of 
their structural parameters. 

In order to circumvent this problem the density pro- 
file was built by superposing RDPs constructed by counting 
stars inside successive radial rings of fixed width up to 55 
arcmin and then dividing by the area of the rings. Ring 
widths ranged from 0.75 — 2.00 arcmin, with 0.25 arcmin 
increments. The narrower rings are ideal to probe the core 
structure of the cluster while the larger ones are better to 
probe the external regions. This prevents either region of 
the cluster to be undersampled. 

The possible contamination from NGC 1977 stars to 
the cluster population was investigated with separate radial 
profiles for the northern and southern region of NGC 1981 
using semi-circles based on the central declination {6 = 
— 04°20^35^^) determined in Sect. 3. Density maps were also 
built up to the extreme radius of 55 arcmin to allow for the 
discrimination of subtle stellar density differences between 
both regions. To obtain a good contrast on the resulting 
density map an optimal radial bin of 2.5 arcmin was used to 
calculate the stellar density around each star. Fig. [5] shows 
the radial profiles of the north (top) and south (bottom) 
regions and the corresponding density maps (insets). 

Although NGC 1977 is about 30 arcmin southwards it 
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R (arcmin) 



Figure 5. RDP of the north (top) and south (bottom) regions 
of NGC 1981. Stellar density maps of each region are also shown 
(insets). Radial ring sizes range from 0.75—2.00 arcmin. Colour- 
bar on the right represents the stellar density (stars. arcmin"^) 
on the insets. 



can be seen in Fig. [5] that the contribution from its stars 
clearly affects the density profile of the southern region of 
NGC 1981. Since it is not possible to tell the boundaries be- 
tween the two clusters or readily distinguish their stars, we 
turned to the northern part of the density profile to inves- 
tigate the structural properties of NGC 1981 by means of 
analytical density functions. 

Trials of King profile fittings were made on the density 
profile using the tw o-parameter modified density function 
introduced bv lKind (^1962): 



/2(i?) = abg + 



ctq 



1 + [R/Rc 



(5) 



To provide better convergence the sky level (crbg) was 
estimated as the mean of the radial bins inside the range 16- 
32 arcmin and further subtracted from each radial ring to 
derive the central density (ctq) and core radius {Re) through 
the 2-parameter King function. The fitting was weighted 
by the statistical Poissonian errors of the star counts inside 
each radial ring. Fig. [6] shows the fitting of the 2-parameter 
function and the sky determination (inset) for the northern 
region of the cluster with the 1-cr uncertainties of the fittings. 
The determined parameters Re and abg are also indicated 
with their associated errors. 

The determined structural parameters were cen- 
tral density cro = 2.4itl.2 stars. arcmin"^, core ra- 
dius Re = 0.83 =b 0.34 arcmin and background level 
Cbg = 0.54 zb 0.07 stars. arcmin"^. A cluster limiting radius 
(^lim) was estimated by visually inspecting the radius 
where stars from the cluster are completely merged with 
the background. We found Rum — 11 ± 1 arcmin for 
NGC 1981. Using the distance of 380 pc (see Sect. 6), 1 
arcmin corresponds to 0.11 pc. 
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Figure 6. 2-parameter King-profile fittings. Background level 
was determined beyond 16 arcmin (inset). The determined core 
radius is indicated. Sky fiuctuation is represented by a grey rect- 
angle. Error bars denote 1-cr Poissonian fiuctuations. 



5 COLOUR-MAGNITUDE DIAGRAMS AND 
FIELD-STAR DECONTAMINATION 

Colour-magnitude diagrams of open clusters usually present 
well-defined stellar sequences such as the main sequence, 
turn-off and giant branch that provide essential informa- 
tion on the physical properties of these objects. However, 
the characterization of open clusters is often hindered by 
strong reddening and field-star contamination, specially for 
the objects projected against the Galactic centre or sparsely 
populated. 

The identification of the cluster stars as opposed to the 
field-star contamination is very important in the study of 
these objects. Assessing membership using proper motion 
data is only available for a small subset of the known open 
clusters, particularly the closest ones. An alternative method 
for disentangling cluster and field-stars consists of using pho- 
tometric data by means of statistical comparison of star sam- 
ples taken from the cluster region and from an offset-field. 
Our decontamination procedure is one of these methods and 
was based on the work of .Bonatto fc B ica (2007. V 

5.1 Decontamination procedure 

A field-star decontamination procedure was carried out by 
firstly selecting an appropriate comparison field presenting 
stellar density and reddening similar to the cluster. While 
circular annular fields circumscribed to the cluster area are 
appropriated for clusters in dense regions and homogeneous 
backgrounds, they fail on regions affected by differential red- 
dening or rapidly changing backgrounds. Therefore, an adja- 
cent field with similar reddening was selected for NGC 1981 
considering the nebulosity contamination from NGC 1977 in 
the southern region of the cluster. Regarding the cluster 
area, the data included stars within the inner 16 arcmin of 
the determined centre. This central region corresponds to 
the cluster visual radius, chosen to encompass the brightest 
B type stars in the field. 
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We built 3D colour-magnitude diagrams for both clus- 
ter and field-stars with J, J — H and J — Kg as axes. These 
colours provide better discrimi nation among the clusters se- 
quences on the CMD (Bonat to et al.ll2Q04l ). The diagrams 
were divided into small cells of average sizes A J = 0.6 mag, 
A( J — H) = 0.3 mag and A( J — K) = 0.3 mag, correspond- 
ing roughly to ten times the average uncertainties in the 
colours. These cells are small enough to detect local vari- 
ations of field-star contamination on the various sequences 
in the CMD, but large enough to accommodate a significant 
number of stars. The grids in Fig. [8] illustrate these cell sizes. 

Initial cluster membership was assigned to clus- 
ter stars within each cell based on their overdensity 
with relation to the field-stars, according to the relation 
P = (A^ciu — Niid)/Nciu- Null probability was assigned when- 
ever an excess of field-stars over cluster stars occurred in a 
given cell. A subset of the original cluster sample was cre- 
ated by removing from each cell on the cluster CMD, the 
expected number of field-stars as measured in the control 
field CMD based on their distance to the calculated centre 
of the cluster. Particularly, each cell on the cluster CMD had 
the stars most distant from the cluster centre removed and 
cells without cluster overdensity had all stars inside their 
limits removed. 

In order to account for the initial choice of parameters, 
we applied the described method for different grid specifica- 
tions by changing the position and size of the cells in each of 
the CMD axes. Cell positions were changed by shifting the 
entire grid one third of the cell size in each direction. Cell 
sizes were increased and decreased by one third of the aver- 
age sizes in each of the CMD axes. Considering all possible 
configurations, 729 different grid sets were used to derive fi- 
nal membership probabilities by taking the average, median 
and mode of the membership obtained in each configuration 
for each star. An exclusion index was also created by noting 
how many times each star was removed from the sample, 
and then normalizing by the number of grid configurations. 

Possible effects of the offset-field selection were inves- 
tigated by applying the method to 3 additional compari- 
son fields 1° distant from the cluster in different directions 
northward of the cluster. The average, median and mode 
of the membership of each star, obtained from the multi- 
ple grid configurations, were compared once different offset- 
fields were employed in the decontamination method. The 
offset-field selection implied mean deviations of 7 per cent 
on the average membership, 8 per cent on the median mem- 
bership and 12 per cent on the mode membership. These 
values were adopted as the general uncertainty of these sta- 
tistical indicators. Final membership probabilities were thus 
assigned to each star by taking the average of the member- 
ship from each grid configuration. 

The decontaminated sample was obtained by removing 
stars based on another two complementary criteria. First, 
stars with average membership value lower than 7 per cent 
or median membership value lower than 8 per cent were re- 
moved. Second, stars with exclusion index larger than 50 per 
cent were removed. The first criterion ensures the removal 
of stars that could present null membership, based on the 
uncertainties of each statistical indicator. The second crite- 
rion is based on the exclusion index and decontaminate the 
cluster sequences by removing the expected field population. 
The exclusion index threshold (50 per cent) was set to pro- 
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Figure 7. Decontamination domain of the complementary 
adopted criteria. Note that stars with lower membership proba- 
bility (left) were removed by the first criterion whereas stars with 
higher exclusion index (right) were removed by the second cri- 
terion. The colourbar indicates both the membership probability 
(left panel) and the exclusion index (right panel) of stars. 



vide the removal of approximately the same number of stars 
removed by the first criterion. 

The adopted criteria complement each other well, act- 
ing on different parts of the CMD and presenting very 
little overlapping of the removed stars. The first crite- 
rion removes stars with low membership probability, acting 
mainly on regions where the cluster and field population 
are entangled, presenting roughly the same density on the 
CMD. For NGC 1981 it represents the lower and mid-main 
sequence/pre-main sequence. For the second criterion, it is 
worth recalling that stars where removed from each cluster 
CMD cell according to the counts from the offset-field CMD 
and on the distance from the cluster centre. Therefore the 
second criterion clean the cluster sequences from the stars 
in the outer regions of the cluster and effectively removes 
from the CMD the stars were field population is dominant. 
It decontaminates the region of the pre-main sequence stars, 
mainly for the fainter regions (J > 10). Fig.[7lshows the de- 
contamination domain of each criterion. 



5.2 Results on NGC 1981 

In order to increase the statistical sample and better probe 
the field population, the procedure was applied on stars 
within 32 arcmin from the cluster centre and the result- 
ing decontaminated sample, used in the subsequent analy- 
sis, was truncated at the visual radius of the cluster (16 ar- 
cmin). This truncation was done because the adjacent south- 
ern clouds associated to NGC 1977 yield a very reddened 
population which occupy a well defined region in the CMD 
(see Sect. 8). The CMD locus of this reddened population 
does not overimpose that of most cluster stars therefore not 
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Figure 8. Cluster (left), field (middle) and decontaminated (right) CMDs using J — H (top) and J — Kg (bottom) colours. Very reddened 
stars, present in the southern region of the cluster were marked with crosses. Stars inside the i^um are marked with dots. The average cell 
sizes are indicated by dotted grids. MS and PMS isochrone fittings are also shown on the right panels. Colourbar on the right indicate 
the calculated cluster membership probability. 



interfering with the decontamination procedure. However, 
even inside the visual radius, very reddened stars present- 
ing J — H > 0.9 and J — Kg > 1.2 appear on the southern 
region of the cluster (crosses in Fig.[8|). They were not prop- 
erly sampled by the field population and therefore could 
neither be excluded by the procedure nor present reliable 
membership probabilities. These stars are likely members 
of NGC 1977. Fig. [9] shows the sky chart of NGC 1981 up 
to 16 arcmin, comparing the member and field population. 
Structural parameters derived in Sect. 4 are also indicated. 

The number of stars removed from the contaminated 
sample do not account for the entire field population be- 
cause the expected number of field-stars measured on the 
offset-field is still higher than the number of excluded stars. 
On average, each criterion removed 45 per cent of the field 
population while presenting a common exclusion rate of 5 
per cent and a total decontamination efficiency of 85 per 
cent. Most of the exclusion takes place in the crowded areas 
of the CMD as 95 per cent of the stars excluded belong to 
the lower main-sequence (J > 12). 

The choice of the offset-field accounts for an average 
deviation of 5 per cent in the number of stars left in the 
decontaminated subsamples (member stars). Moreover the 
fraction of recurring member stars, independent of the offset- 
field, is greater than 87 per cent. Fig. [8] compares cluster, 
field and decontaminated stars in the J x J — H and J x 
J — Ks diagrams. 
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Figure 9. Sky chart comparing the removed field population 
(open circles) and remaining cluster stars (filled circles). Concen- 
tric circles represent the core radius and limiting radius respec- 
tively. Very reddened stars are marked with crosses. Colourbar 
indicates assigned membership probabilities. 
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Figure 10. Correlation between our membership assignment and 
proper motion data taken from UCAC2 (top). Mean proper mo- 
tion uncertainties are indicated. Probability distribution function 
of member and non-member stars are also shown in the bottom 
panels for RA (left) and DEC (right) proper motion directions. 

5.3 Comparison with proper motion data 

Proper moti on data from UCAC2 ([Zacharias et al.l |200J) 
and UCAC3 (Zacharias e t al.ll2010l l catalogues were used to 
investigate the reliability of the photometric memberships 
determined for the cluster stars. We selected stars within 
15 arcmin from the cluster centre, with 2MASS photom- 
etry subject to the constraints described in Sect. 2.1 and 
only those with proper motions and their uncertainties. The 
corresponding samples amount to 159 objects in UCAC2 
catalogue and 99 objects in UCAC3 catalogue. 

Comparison of the UCAC2 proper motion with our 
membership results was done using the 2MASS designation 
to identify the common stars between the datasets. This al- 
lowed us to use our photometric membership to discriminate 
cluster and field-stars on the vector point diagram (VPD) 
(Fig. 1101 top panel). Although the photometric members 
present a concentration near (|j,RA, |j,DEC) = (0,0), this is 
also true for the field stars, as can be more easily seen in the 
probability distribution histograms (bottom panels). 

As the statistical distance between the cluster and field 
proper motion distributions is small, it is difficult to sep- 
arate the populations. This effect becomes more accentu- 
ated for clusters where there is a large field to cluster 
member ratio or when the centroids of the proper motion 
distributions are very close to each other, such as in the 
case of NGC 1981. This problem has already been discussed 
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Figure 11. Correlation between the proper motion membership 
from DALM02 and our photometric memberships (left). Stars be- 
tween the dotted lines have membership differences smaller than 
20 per cent (right). 



by Cabrera- Cafio fc Alfarol (Il990) and Sanchez fc Alfarol 
(2009) which proposed to use the spatial distribution of the 
stars, in addition to their proper motion distribuition, to 
increase the statistical distance between the populations. 

In fact, membe rship probabiliti es have been estimated 
for NGC 1981 by Bias et al. (2006) using UCAC2 proper 
motion of 160 stars within 15 arcmin of the cluster centre. 
The comparison between these memberships and the photo- 
metric memberships derived in this work shows a very poor 
correlation, as can be seen in Fig. [TT] (left panel). Although 
the small statistical distance between the populations may 
account for the large scattering in the proper motion mem- 
bership, especially for the low-membership stars, it seems 
that the derived photometric memberships are, in average, 
underestimated in relation to the proper motion ones. This 
trend is still visible when the low-membership stars are re- 
moved from both the photometric and proper motion sam- 
ple, although the correspondence between the samples im- 
proves as most of the remaining stars present membership 
differences smaller than 20 per cent (right panel). 

UCAC3 provided a complete new reduction of the 
Southern Proper Motion data that improved the proper mo- 
tions of faint stars by a factor of 2 compared to UCAC2. Al- 
though this first releas e still presents some unsolved prob- 
lems as described in Zacharias et al.l (|2010|), we present a 
comparison of our results with proper motion from UCAC3 
in Fig. 1121 It is clear that UCAC3 VPD is sparser and less 
populated than that using UCAC2 data. However, it is cur- 
rently not po ssible to realize the origin of this lack of stars 
(see Sect. 8 in lZacharias et al]|2010 V Additionally, while our 
photometric analysis yields very entangled cluster and field 
populations in UCAC2 VPD, they are much better discrim- 
inated when proper motions from UCAC3 are considered 
(bottom panels of Fig. [TOl and [T2|) . Therefore, membership 
derived from UCAC3 proper motions might be in closer 
agreement with our photometric analysis. 



6 ASTROPHYSICAL PARAMETERS 

The optical data was correlated with the infrared data al- 
lowing the membership probabilities derived from 2MASS 
data to be assigned to the stars with BVRI photometry. 
Although shallower than the infrared sample and presenting 
larger photometric uncertainties, the optical data provided 
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Figure 12. Similar to Fig.[Tol but using UCAC3 proper motions. 



a larger spectral base for the CMD analysis and additional 
constraints on the physical parameters. 

Astrophysical parameters were determi ned by means 
of f ittings of Padova isochrones fo r 2MASS (Bonatt o et al.l 
12004 ) and Johnson- Cousins ( Marigo et al.l i2008 ) filters on 
the dec ontaminated CM Ds. Pre-main sequence evolutionary 
tracks (I Siess et al.ll200 a) were also employed. We chose only 
isochrones with overshooting and solar metallicity. As initial 
values for the cluster physical para meters we wer e gui ded by 
the pr evious studies of .Kharchenko et al.l (|2005|) and iBallvl 

liooi). 



6.1 Age 

The absence of evolved sequences in the CMDs due to the 
cluster youth makes determination of its age unreliable if 
one uses standard evolutionary models, i.e., those with- 
out accounting for pre-main sequence stars. Pre-main se- 
quence isochrone fittings suggest that stars with ages rang- 
ing from 1-10 Myr coexists in the cluster and that most 
high-membership (P > 0.6) stars seems to fall between the 
young pre-main sequence isochrones of 1-5 Myr. This be- 
haviour is also observed when only the innermost stars (in- 
side the limiting radius) are considered. Fig. [8] shows the 
isochrone fittings on the decontaminated J x (J — H) and 
J X (J — Kg) CMDs along with the derived membership of 
stars. Fig. [T3l shows the same fittings on the V x {V — I) and 
V X (V — J) CMDs. The average age derived for the cluster 
is 5± 1 Myr based on isochrone fittings to both near-infrared 
and optical data. 




Figure 13. NGC 1981 decontaminated CMDs with optical data. 
Colours and symbols are the same as in Fig. [S] 



In photometric studies of star clusters the MS turn-off 
is the most reliable feature for determining the age of these 
objects. However the identification of this CMD feature in 
very young clusters is uncertain by the scarcity of massive 
stars and the isochrones lack of turn-off sensitivity. For these 
objects the turn-on, the CMD locus where the PMS stars 
joins the MS, is a good indicator of the age of the stellar 
population. 

From the isochrone fittings shown in Fig. [8] it can be 
seen that NGC 1 981 turn-on is at J ^ 9, corresponding to 
an age of 5 Myr. Cigno ni et al.l ([2010) devised a method to 
derive the age of young stellar systems by detecting the PMS 
turn-on with the luminosity function (LF) . The clustering of 
stars at magnitudes near the turn-on produces a bump in the 
LF that becomes increasingly fainter as the cluster ages. By 
using synthetic simple stellar populations they calibrated a 
relation between the magnitude My of this luminosity bump 
and the age of the population. 

The reddening and distance modulus determined in the 
present work were used to construct LFs with Mj and My 
for both the original and decontaminated sample, which are 
shown in Fig. [M] In both cases it was possible to identify a 
smah bump in the LF at Mj = 0.8 ±0.2 and My = 1.3 ±0.6. 
This bump on the LF correspond to the age 3.3 ± 2.7 Myr, 
according to the age calibration with My bv lCignoni et al.l 
(|2010i). ^__^_^ 

However, as stated by ICignoni et al.l (|2010|), reliable 
identification of the turn-on can only be made in populous 
clusters containing at least ~ 50 stars brighter than My = 5 
in order to successfully detect the bump among the Poisso- 
nian fluctuations. Therefore, although we should use this 
result with caution, it is in agreement, within the uncertain- 
ties, with the age and distance modulus obtained by means 
of isochrone fittings. 
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Figure 14. Luminosity functions in the infrared (top) and optical 
(bottom) bands. The turn-on bumps at My =1.3 and Mj =0.8 can 
be identified on both original (left) and decontaminated (right) 
samples. 



6.2 Distance and Reddening 

Although there were small discrepancies in the de- 
rived parameters (mainly on reddening) from optical 
versus infrared data, the mean determined values of 
distance modulus J — Mj = 8.0 ± 0.1 and colour excess 
E(B — V) = 0.07 =b 0.03 provided good isochrone fitting in 
both cases. With those values the true distance modulus is 
(m - M)o = 7.9 ± 0.1 whe re we have used Aj/Ay = 0.282 
and Ay/E{B -V)^ 3.09 (JRieke k Lebofskvlll985l ). result- 
ing in the distance of 380 =b 17 pc. The uncertainties were 
set to accommodate the values derived from the optimal fit- 
tings on both optical and infrared CM Ds. Thes e values are 
similar to the values found by Kharche nko et al.l 12005). The 
almost 50 per cent deviation between optical and infrared 
reddening values deserves further discussion. 

Although the literature reddening value provided a 
good isochrone fitting on the infrared data, the reddening 
value obtained from the V x {V — I) CMD suggests that the 
extinction in the optical bands is larger by a factor of 2. 

Reddening in the o ptical bands was also determined in 
a previous w ork (Maial boOTl l by a hnear fit to the ZAMS 
(ISchmidt-Kalerll 198211 in co lour- colour diagrams as proposed 
bv lMunari &: Carrarol (|l996, V By selecting the diagram most 
sensitive to reddening (Fig. [15)) and limiting our analysis to 
stars brighter than y = 13 we derived a reddening value of 
^(^ — y) = 0.11, which is very close to the value obtained 
from our best isochrone fitting on the optical data. 

The reddening i n the cluster di r ection , as interpolated 
in the dust maps bv ISchlegel et al.1 (|l998l \ is E{B - V) = 
0.24 ± 0.04. Since these maps are sensitive to dust column 
density, this value can be considered as an upper limit to the 
cluster reddening. Also the cluster may be affected by dif- 
ferential reddening caused by the dust clouds in the region, 
specially for its southern part. 



7 RADIAL MASS FUNCTION 

To analyse the spatial dependence of the mass function 
(MF), the decontaminated LF in the J band determined 




Figure 15. Colour-colour diagram with the optical data showing 
the ZAMS. 



in Sect. 6.2 was used to investigate the regions compre- 
hended within the inner 5.5 arcmin (half limiting radius) 
and within the range between 5.5 arcmin and the limiting 
radius (11 arcmin) of NGC1981. The decontaminated LFs 
were converted into one MF for each region by fitting a mass- 
luminosity relation (in J band) from the combined PMS and 
evolved isochrones of 5 Myr at solar metallicity. The power 
law (/)(?7i) = A • m~'''^~^^\ where A is a normalization fac- 
tor and X is the MF slope, was fitted to the data in both 
regions. The MFs were normalized by the number of stars 
(N) inside each region (N — 25 for the inner region and 
N = 81 for the outer region). The outer region MF fit was 
limited to stars less massive than ?^ 1 M© due to its small 
number of brighter stars. The MFs of the selected regions 
are presented in Fig. 1161 together with the power-law fits 
and the determined slopes. 

Both regions have MFs that are flatter than the lSalpeterl 

(|l955h one, with slopes x — 0-65 =b 0.08 (inner region) and 
X = —0.44 zb 0.03 (outer region). Total mass inside each 
region was calculated by summing over selected mass ranges. 
We obtained for the inner region m = 30 ± 6 M© (0.4 < 
m(M0) < 6) and for the outer region m — 107 ± 13 M© 
(0.2 < 771 (M©) < 9). For the inner region the mass was 
also calculated by integrating the derived power law over 
the same mass range yielding m = 38 =b 2 M© . We should 
note that this mass implies an uniform distribution of stellar 
masses while the observed distribution of the brightest stars 
suffers from small number statistics. 

The relaxation time of a star system can be defined as 



^relax 



ytc 



s, where tcross = R/ctv is the crossing time. 



N is th e total number of stars and ay is the velocity dis- 
persion teinne v fc Tremaind [l987l : iLada fc Ladal 120031 ). By 
using the dispersion in proper motion given by DAML02 for 
31 stars within 11 arcmin from the cluster centre with mem- 
bership probability (derived by DAML02) above 50 per cent 
and assuming an isotropic distribution of spatial velocities 
we estimated a mean velocity dispersion cr^ ~ 8 km.s~^. 
Combining this information with the limiting radius we get 
tcross = R/(Tv = 0.15 Myr and treiax = 0.41 Myr for N = 106 
stars. 

Although unexpected for such a young object, dynami- 
cal evolution can be responsible for rapid collapse of the clus- 
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Figure 16. Mass function of stars inside the inner 5.5 arcmin 
(top) and within the range 5.5-11 arcmin (bottom). Power law 
fittings, resulting slopes and total mass are indicated. Error bars 
correspond to 1-cr fluctuations 



ter core causing mass segregation and evaporation of stars 
in outer regions. This may have occurred in NGC 1981 since 
the relaxation time of the cluster is much smaller than its 
age. As shown in Fig. [16] the mass function slope is larger 
for the inner 5.5 arcmin than for the outer 5.5-11 arcmin 
annular region between the mass range 0.4-1.1 M©. This is 
in apparent contradiction with the expected consequences 
of dynamical evolution, i.e., the existence of an excess of 
low mass stars and depletion of massive ones (eventually 
grouped in binaries and/or multiple systems) in the outer 
region compared to the inner region. However we indeed do 
not found any star with mass lower than 0.4 M© in the in- 
ner region but 23 stars in the same mass range in the outer 
region, in agreement with the mass segregation scenario. 



8 STAR DENSITY MAPS 

We have mapped the projected spatial star density and the 
frequency of stars in different positions of the CMD (Hess 
diagram). Firstly we show in Fig. [TJlthe star density map 
and the corresponding Hess diagram of a large region around 
NGC 1981 including the Orion Nebula Cluster (ONC). Both 
diagrams were built by dividing the plots in square cells of 32 
pixels on a side, which in Fig. [TTl corresponds to 5.6 arcmin 
for the spatial map, 0.1 in (J — Kg) and 0.4 mag in J for 
the Hess diagram. The density diagrams are then smoothed 
by a cubic convolution interpolation method. 

The data comprises all sources extracted from 2MASS 
catalogue obeying the photometric quality constraints men- 
tioned in Sect. 2.1. The ONC dominates the density map 
southwards of NGC 1981 which the centre, as determined in 
Sect. 3, is indicated by a square. Even though, an overden- 
sity compared to the adjacent neighbourhood can be noticed 
at the position of NGC 1981. The Hess diagram provide ev- 
idence that most of the stellar content in the ONC consists 
of red low mass stars. 

Smaller sky and CMD areas were considered around 
NGC 1981 in Fig. [18] The square cells of 32 pixels correspond 
to 2.4 arcmin for the spatial map and 0.02 in (J — Ks)^ 
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Figure 17. Star density map (left) and Hess diagram (right) of 
a region 3°x3° centred in NGC 1981 (square) based on 2MASS 
data. 
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Figure 18. Same as Fig. [TT] but for a smaller region centred in 
NGC 1981 and (J - Ks)< 0.6. 



0.3 mag in J for the Hess diagram. The spatial map has 
1.3° X 1.3° and the associated Hess diagram has been limited 
to contain only stars with {J — Kg) < 0.6, thus excluding the 
reddest stars. Such a colour constraint clearly separates the 
spatial map in two regions: a denser star field northwards of 
NGC 1981 (but including it) and a more rarefied star field 
southwards of NGC 1981. It clearly reflects the presence of 
dust revealing a boundary between a heavy obscured stellar 
field and a region with lower extinction. It also indicates 
where the reddest stellar population lies in the area. The 
cut in colour eliminates reddened stars by dust but also pre- 
main sequence stars, known to be abundant in the ONC 
region. 

Fig. [19] shows the spatial stellar density plots for the 
same area as in Fig. [18] but with progressive cuts in star 
brightness, allowing one to connect the density peaks with 
different populations characterized by a limit magnitude. 
The plots show the emergence of NGC 1981 as the limiting 
magnitude varies from J= 15.8 to 10. 

Although NGC 1981 is in a region of intense star forma- 
tion with clouds of dust and gas, it has presently inhabited 
a relatively dust-free field, perhaps as a consequence of the 
cluster evolution. Its massive stellar content may have con- 
tributed to the energy release into the interstellar medium, 
either by means of exploding supernovae and/or winds due 
to radiation pressure. Our study concludes that NGC 1981 
is a cluster older than the ONC and NGC 1977, which are 
still embedded in the parental gas and dust cloud, and that 
its structure is best visualized by considering stars brighter 
than J = 12. 

Fig. [20] shows the same analysis applied to probable 
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Figure 19. Star density maps for (J — Ks)< 0.6 and variable J 
magnitude limit as indicated. 
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Figure 20. Decontaminated star density map (left) and Hess 
diagram (right) for a circular region oi R = 15 arcmin centred in 
NGC1981. 



members inside a circle of radius i? = 15 arcmin centred in 
the cluster. As a consequence of the method applied for field 
decontamination, the cluster stellar content is more clearly 
defined down to the photometric limit of the data. The un- 
derlying stellar field and the cluster stellar population were 
indistinguishable in the previous contaminated diagrams for 
stars fainter than J = 12 (Figs. [n]to[19l). In Fig. [20l a 
population of pre-main sequence members stands out in the 
Hess diagram. 



9 DISCUSSION AND CONCLUSIONS 

We investigated the stellar object NGC1981, a young star 
cluster in the Gould's Belt nearby the Orion Nebula Cluster 
containing a small population of massive B type stars. We 
employed near-infrared data from 2MASS to calculate the 
cluster centre at a = 5^35^08", 6 = -04°20'35'' and to 
subsequently derive the structural parameters from King- 
profile fittings obtaining a core radius Re = 0.09 dz 0.04 pc 
and central density ctq = (2±1) x 10^ stars. pc~^ . A limiting 
radius of Rum = 1-21 =b 0.11 pc was also derived from the 
radial density profile. 

We devise d a decon tamination procedure based on the 
method by Bo natto fc Bica (2007) to statistically remove 
the underlying field population from the cluster CMD by 
using an offset-field to sample the background contamina- 
tion. Tested on NGC 1981 with multiple offset- fields, the 
procedure reliably presented an average field-star exclusion 
efficiency of 84 per cent, separating cluster members with 
an average deviation of 5 per cent in the number of stars. 
It also provides photometric membership for the member 
stars with average uncertainty of 6 per cent, depending on 
the selection of the offset-field. 

By using optical BV{RI)c data alongside 2MASS data 
we performed isochrone fittings on the decontaminated data 
and determined average values of reddening E{B — V) — 
0.07 =b 0.03, distance modulus {m - M) ^ 7.9 ± 0.1 {d = 
380 ±17 pc) and age 5 =b 1 Myr. The rela tion between the 
turn-on My and age (jCignoni et al.ll2010l l was also used as 
an additional check for the cluster age. In this case the age 
derived was 3.3 ±2.7 Myr. The scarcity of bright stars make 
this value unreliable as sole indicator of age. We use it to 
infer that the luminosity of the turn-on is consistent with 
the age and distance modulus found by isochrone fittings. 

By comparing our derived memberships with proper 
motions from UCAC2 catalogue, we showed that while the 
small statistical distance between cluster and field popu- 
lations may have hardened their separation in the VPD, 
the photometric method was capable to discern these very 
entangled populations. Indeed, by comparing our derived 
memberships with proper motions from UCAG3, we were 
able to distinguish a clear concentration of member stars 
over a sparse distribution of non- member stars in the VPD. 

We derived mass functions for stars inside the inner 
5.5 arcmin and 11 arcmin (limiting radius) and evaluated 
the total mass inside these regions: ?7i = 30 ± 6 M© and 
m = 107 ± 13 M0, respectively. After fitting a power law 
to the data, the calculated slopes x — 0.65 ± 0.08 (inner re- 
gion) and X — — 0.44±0.03 (outer region) are flatter than the 
Salpeter MF, indicating a depletion of low mass stars in the 
outer region relative to the inner region. Cluster evolution 
and mass segregation can explain this effect if the cluster has 
had time to dynamically evolve. In fact, by using proper mo- 
tion data from DALM02 for stars inside the limiting radius 
we have derived tcross = 0.15 Myr and treiax = 0.41 Myr, 
meaning that the cluster stars had time to interact gravi- 
tationally with each other forming binaries and/or multiple 
systems composed by massive stars that settle in the cluster 
centre and leading to the ejection of low mass stars towards 
the cluster outer regions. 

The cluster might have undergone mass loss in the ear- 
lier stages of the star forming process leading to the present 
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sparse structure. The energy released into the medium by 
the evolution of the most massive stars (i.e. by supernovae 
and stellar winds) may have cleared the cluster from its 
parental cloud of gas and dust causing a collapse of the core 
due to the changing gravitational potential and subsequent 
evaporation of stars in its outer regions. 

Through stellar density maps we discriminated the pres- 
ence of a embedded red population of stars just south of 
NGC 1981. These stars are likely members of the young ob- 
ject NGC 1977 and still contaminate the southern popula- 
tion of NGC 1981. These maps also demonstrate the power 
of the procedure employed to disentangle field-stars from the 
cluster population. 

As the majority of young open clusters, NGC 1981 is 
not expected to live longer than a few Myr, evolving from 
its actual state of marginally bound system to a loose stel- 
lar association and finally dispersing itself into the Galactic 
disc. 
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